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We study Bs -^ KK decays in the framework of the PQCD. We show that the branching 

^ '■ ratios of 5, -^ {K+K',K°K^) are about (22.43,25.78) x lO^^ for ^3(7) ^ 72°, which are 

consistent with the model-independent estimations. We find that the typical scale of the 

decays is near 1.7 GeV. We also point out that the induced strong phase 6 is around 207° 



so that the direct CP asymmetry of Bs — > K^K could reach 15% 



0. 



The study of charmless B decays has an enormous progress since many decays such as 
those with exclusive pseudoscalar-pseudoscalar final states {B — > PP) were measured at 
e~^e~ machines by CLEO [jl| , B AB AR p| and Belle 0] , respectively. From the search of B 
decays, we not only could test the origin of CP violation in standard model (SM), which is the 
consequence of the Cabibbo-Kobayashi-Maskawa (CKM ) quark-mixing matrix 0] but also 
verify various QCD approaches for nonperturbative problems in exclusive decays. Recently, 
the proposals of using the mixing-induced and direct CP asymmetries in B^i -^ tt+tt^ and 
Bg -^ K^K^ [^ or the branching ratios (BRs) of Bg —>■ K^K~ and K^K^ decays [^ are 
suggested to determine the angle 03 or 7. Clearly, it is important to give a detailed analysis 
on the decays of Bg — > KK. 

It is known that one of the main theoretical uncertainties for studying exclusive hadron 
decays is from the calculations of matrix elements. Usually it is performed in the pertur- 
bative QCD (PQCD) approach developed by Brodsky and Lepage (BL) |^. In the BL 
formalism, the nonperturbative part is included in the hadron wave functions and the tran- 
sition amplitude is factorized into the convolution of hadron wave functions and the hard 
amplitude of the valence quarks. However, with the BL approach, it has been pointed out 
that perturbative evaluation of the pion form factor suffers a non-perturbative enhancement 
from the end-point region with a momentum fraction a; ^ p[. If so, the hard amplitude is 
characterized by a low scale, such that the expansion in terms of a large coupling constant as 
is not reliable. Furthermore, more serious end-point (logarithmic) singularities are observed 
in the twist-2 (leading- twist) contribution to the B ^ ir transition form factors ^, p!o[| . 
The singularities even become linear while including the twist-3 (next-to-leading twist) wave 
function ||rT| . Because of these singularities, it was claimed that form factors are dominated 



by the soft dynamics and not calculable in the PQCD [|T^. 

In order to take care of the end-point singularities, /c^, the transverse momentum of the 
valence-quark []T3|, and threshold resummations []T1|, |I^] have to be introduced. The inclusion 



of kx will bring in large double logarithms as\yi^{kT/MB) through radiative corrections. 
Therefore, these large logarithms should be resumed in order to improve the perturbative 



calculation. Due to the resummation |jT3|, |T6|, |T^, the arisen distribution of kx exhibits a 
suppression in the region with k^ ~ 0(A) ||19| and the average of k^ is enlarged up to around 
(k^) ~ AMb for Mb ~ 5 GeV; consequently, the off-shellness of internal particles keeps 
being 0{AMb) even in the end-point region. Moreover, due to the radiative corrections of 
the weak vertex |T^, another type of double logarithms ag In x actually exists while x ^ 0; 
and therefore, these large corrections should be also resumed, called threshold resummation 



I^, O, for justifying the perturbative expansion so that the result leads to a strong Sudakov 



suppression at x — > 0. Hence, including kx and threshold resummations, the end-point 
singularities can be dealt with self-consistent in the PQCD. 

In the literature, the applications of the modified PQCD (MPQCD) p2[ to the processes 
of fi ^ PP, such as B ^ K-n [|T1, fi ^ tttt pl, fi ^ KK H and 5 ^ K^'^ [||], 



as well as that oi B -^VP such as 5 ^ (/)7r [gg, 5 ^ p(cj)7r [g^ and 5 ^ (t)K @, gg 
have been studied and found that they are consistent with the experimental data or limits. 
For a review on the PQCD approach, we summary the characters of the MPQCD briefiy as 
follows: 

• Due to the introduction of kr and threshold resummations for smearing the singular- 
ities, the B ^ 7T,K form factors are still dominated in the perturbative region with 
a,/7r<0.2 pq. 



There involve three scales in the MPQCD: the My/ scale of the electroweak (EW) 
interaction, the typical scale t which reflects the speciflc dynamics of the heavy meson 
decays, and the factorization scale of A ~ Mb —Mb with Mb and M^ being the B-meson 
and b-quark masses, which shows the interface of the soft and hard QCD dynamics, 
respectively. We note that the t scale is chosen such that the contributions from higher 
order effects are as small as possible [^ and one can flnd the speciflc scale to be the 
chiral symmetry breaking scale pll, E9|. 



Penguin enhancement: As known, Wilson coefficients (WCs) oiC^i^jji) and Cq^ii) gen- 
erated from the QCD penguin [^ increase significantly at t < MbI'2.. Due to the 
enhancements, it is realized that the BRs of -B — > (\)K in the MPQCD |2^, |29 
explain the results of Belle pOf and BaBar [Rll] naturally. 



can 



Less theoretical uncertainties: The main theoretical uncertainties in the MPQCD from 
are the shape parameter lob of the B-meson wave function, chiral symmetrical breaking 
parameter Wip introduced by the matrix element of the pseudoscalar nonlocal oper- 
ator for the P meson |T^, ^, and the power factor c for the parametrization of the 
Sudakov factor St oc [a;(l — x^f generated by the threshold resummation, respectively. 
With experiments, one would precisely set the limits on these theoretical unknown 
parameters. 

Large absorptive parts: The major absorptive parts arise from the annihilation topol- 
ogy in which the propagator of the internal quark satisfies the on-shell condition. With 
the power counting rule, the ratios of the transition form factor (F^^) to the imaginary 
and real parts of annihilation contributions are found to be F^^: ImF^^: ReF^^ 
= 1 : 2m%/MB: A/ Mb (H- For Mb ~ 5.0 GeV, the value of the imaginary part is 



compatible with that of the form factor. 



In this paper, to calculate the matrix elements of four-quark operators relevant to the 
Bs -^ KK decays, we adopt the MPQCD factorization formalism as 



[dx\ 



An 



^*k{xM^k{x2M) 



{VK\Ok\B) 



X 5,({a:})e~"V— --— y (1) 

where $^ and $Bs ai's the wave functions of K and Eg., T^ is the hard scattering amplitude 



X Tk{{x],{h],MB)^BAxiM 
-5({x},{6},Mb 



dictated by Ok operators, the exponential factor is the Sudakov factor |jT6|, |T7|, and St{x) 
p!8| , p!9|] expresses the threshold resummation factor. Since B^ is a heavy meson, the chiral 
symmetry breaking effects are negligible so that ^b^ is regarded as $Bj„. 
The effective Hamiltonian for decays with the b —* s transition is given by 



F, 



G 



cff 






Ci{fi)o[' ^ + C2(/i)o^^ ^ + E Ci(/^)^ 



(g') 



10 



j=3 



(2) 



where Vq' 



VgigVq'b are the products of the CKM matrix elements, Cj(/x) are the WCs 
and Oi correspond to four-quark operators, and their detailed expressions can be found in 



Ref . P2| . It is easy to find that the structures of 03,4,5,6 , generated by the QCD penguin, are 



the same as that of 09,10,7,8 from the EW penguin, respectively. Except the different color 
flows between 02i-i and O2J {i = 1, 2, ■ ■ ■ , 5), each of the pair operators also has the same 
structure. Therefore, we can define the useful EW dynamical variables as 

_ r^ , ^2 _ r^ , ^1 ' _ ^2,1 

ai - ^i + J^^^ ^2-62 + —, «i,2-^, 

%,4 — <-'3,4 + -^<-^9,10 + «3,4; «3,4 — "^ + ^^^-'lO.Q , 

where the superscript q can be the light u, d and s quarks, respectively. We note that all 
decay amplitudes related to the weak dynamical interactions are dependent on of and a^ 
(j = 1, 2, ■ ■ ■ , 6), while nonfactorizable effects are only associated with the color suppressed 
variables d^ 



'(g) 
By defining the decay rates as 



327r ' ' 

with Gp and Mb being the Fermi constant and the mass of Bg, the corresponding decay 
amplitudes for B^ -^ K^K" and Bg —>■ K^K^ modes are written by 

A^- = fKv:F:^^^ + v:MT + fBv:F:,^'^ + v:{MT+MT+MT) 

-fKv:F,, - v:Me2 - v:Mai, 
A'' = fKv:F:^^^ + v:M:i^ + fBv:F:,^'^ + v:{MT + MT+MT), (4) 

respectively, where fK{B) is the K{B) decay constant. In Eq. (|), {F,^^} denote the 
factorizable emission (annihilation) transition form factors for tree [penguin] diagram con- 
tributions, while {A^e((j) } correspond to nonfactorizable effects, the superscript q represents 
the q quark pair emitted from the EW penguins, the subscripts of 1-6 label the WCs of Eq. 
(0) appearing in the factorization formulas, and the nonfactorizable amplitudes M.e2(ai) are 
from the operators Oi^l- Note that there is no tree contribution to the decay of Bg -^ K^K^ . 
Explicitly, one has that 

pP(g) _ pP{g) pP(g) 

F^^''' = SttCfMb / dxidx2 / bidbib2db2(t)Bixi,bi) 
Jo Jo 

X {[(1 + X2)0A'(1 - X2) + rKil - 2X2) (0^(X2) + 0^(a;2))] 
xE^J(t«)/le(Xi,X2,6l,fe2) 

+2rK<p''K{x2)Eif{t'i^)K{x2,x,,b2,b,)} , (5) 

Fg ^^^ = IGnCpM^^rx / dxidx2 / bidbib2db2(l)B{xi,bi) 

Jo Jo 

X {[0/c(l - X2) + Tk ((2 + X2)0^(X2) - X2(I)k{x2))] 

xEif{t'^^^)K{x,,X2,bi,b2) 

+ [2rKfK{x2)]E^\t^^^)K{x2,Xi,b2,b,)} , (6) 



F„g = IGnCpM^^rK / dx2dx^ / h2dh2h^dh^ 

Jo Jo 

xE^j){t'i^)K{x,,X2,h,b2)} , (7) 

with rx = m%/MB for q = u- or rf-quark. Here, Cp = 4/3 is the cofor factor, (pB and 0^ 
befong to twist-2 Bg and K wave functions, while 0^ and 0^ are for twist-3 [Q and their 
detailed expressions can be found in Ref. [^, the hard part functions he(a) have included 
the St factor [^ and the evolution factors are given by 

Ei^\t) = a,{t)a^\t)exp[-SB{t,Xi)~SK{t,X2)], (8) 

Ei^\t) = a,{t)a^\t)exp[-SK{t,X2) - SK{t,Xs)], (9) 

respectively. As expected that the nonfactorizable effects are smaller and more complicated, 
we do not display their expressions here but they will be included in the numerical calcu- 
lations. We note that in our PQCD approach the x dependence in the kaon wave function 
is usually assigned to the momentum fraction of the light u oi d quark ||29| . However, due 



to the s quark being spectator in emission contributions, for simplified the hard parts, we 
adopt that this s quark carries a momentum fraction of x so that in order to compensate 
this change, the arguments X2 of the kaon wave functions {(J)k (3^2)} in Eqs. (||) and (^ are 
replaced by 1 — X2 in which we also have used the properties of 0^ (1 — X2) = 0^ (^2) and 
01^ (1 — X2) = —(pK (^2) 1^, P^- On the other hand, the factorizable annihilation contribu- 
tion, associated with WC a^ and described by 



F„^^^^ = IGnCpMl f dx2dx3 H h2dh2hdh 

Jo Jo 

X { a;30x(a;2)0i^(a;3) + 2r^0^(a;2) ((1 + a;3)0^(x3) 

-(1 - X3)0^(X3))] 4J(t«)/l„(x2, X3, 62, &3) 
- X2(t)K{x2)'pK{x-i) + 2r^0^(x3) ((1 + X2)(pKi^2) 

-(l-X2)0^(x2))]4J(tf)/i„(x3,X2,63,&2)} , (10) 

vanishes. It can be seen easily by interchanging the integration variables X2 and X3 in the 
second terms of Eq. ([T0|) . However, this property does not apply to the annihilation contri- 
butions associated with Og which are constructive in Eq. (^. The factorization formulas 
for Fai and F^^^ , associated with the WCs of ai(ta) and 03 (t^) + 05 (ta), are the same as 
Fa4 , i.e., vanish. 

In the numerical analysis, we first show the Bg ^ K form factor of F^"^ in Figure jT], 
which can be easily obtained by dropping the WC dependence out of Eq. (^, where we have 
used fB = 0.20 and /x = 0.16 GeV. For a comparison, we also display the results oi B^ ^ K 
form factor i^F^'^^) from {(px (3^2)} and Mb^ instead of {0x (1 — X2)} and M^^, respectively. 
From the figure, we see that F^"^ is slightly less than F^-^^ and the main difference is from 
the chiral symmetry breaking effect appearing in (pK arisen from the different argument of 
(pK- After including WCs, the values of individual terms in Eq. (^ can be read from Table ^ 



Table 1: The values of individual transition form factors (TFFs) for Bs -^ KK decays with 



Ub 



0.4, two sets of m?^ and A^^,"^ = A^fi'^ + A^Sb^ + M^^"^ 



'a46 



TFF 


-^646 


-^646 


pP{s 


Fe2 


y\/le46 


m%^ = 1.7 


-31.91 


-32.67 


-1.25 + ^10.91 


369.37 


0.17 + i0.32 


m5^ = 1.5 


-27.00 


-27.64 


-1.10 + ^9.63 


336.53 


0.16 + i0.28 


TFF 




mT 


My 


Me2 


Mai 


m°^ = 1.7 


0.16 + i0.48 


0.19 + i0.52 


0.20 + i0.54 


-1.19-22.75 


0.69 -i4.3 


m?f = 1.5 


0.16 + i0.45 


0.15 + i0.50 


0.15 + i0.52 


-1.36-i2.95 


0.59 - i4.04 



The slight difference between F^/^q and Fg4g conies from the EW effects. The factorizable 
annihilation and nonfactorizable contributions are complex because the on-shell condition 
can only be satisfied in these diagrams. On the other hand, in order to know the strong 
phase more clearly, we reparametrize the decay amplitudes as 



A 



Y*P _ Y*T 

t u 



-v:p 



1 + 



K,. 



Vr 



re 



i{S+<t>3) 



(11) 



with re = 1 + T/P, where we have used Y.i 



-{u,c,t} 



Vi = 0, T and P express the whole tree 



and penguin contributions, 6 describes the strong phase, and the values of r and 6 are shown 
in Table ^ From the table, we clearly see that cos6 < in the MPQCD approach. 



Table 2: The results of r and 6 with different values of ub and m^. 



< 


1.7 


1.6 


1.5 


UJb 


r 


5(deg) 


r 


5(deg) 


r 


5(deg) 


0.41 


9.16 


207.7 


9.42 


208.0 


9.70 


208.3 


0.40 


9.22 


206.9 


9.49 


207.2 


9.78 


207.6 


0.38 


9.34 


205.3 


9.63 


205.6 


9.92 


205.9 



To calculate the CP average BRs of B^ 
values: 



KK decays, we use the following data as input 



K,,. 



A, K, ^ -A\\ K 



ub 



AX^RhC 



-i4>3 



A ^ 0.80, A ^ 0.22, Rh ^ 0.36. 



From these values, the results with possible ub and m^ are displayed in Table |^. We note 
that although the transition form factor of the tree contribution is larger than others by over 
one order of the magnitude as shown in Table [^, due to the suppression of CKM matrix 
elements, actually the tree contribution is subdominant. By fixing WCs at some specific 
scales, the results are given in Figure |^. From that, we clearly see that the typical scale for 
the Bs -^ KK decays is around 1.7 GeV. In Figure ^ we also show the BR of Bs -^ K^K~ 
as a function of angle 03. 

As usual, the direct CP asymmetry (CPA) can be defined as 



A, 



CP 



T{Bs -> K+K-) - T{Bs -^ K-K+) 
T{Bs -^ K+K-) + T{Bs -^ K-K+) ' 



Table 3: The BRs (in units of 10 ^) of Bs —>■ KK decays with different values of ujb and 



m 



K- 



Br 


B, ^ K+K- 


Bs -> K^K^ 


m\ = 1.7 


m% = 1.6 


m^ = 1.5 


m\ = 1.7 


m*^ = 1.6 


m% = 1.5 


ujb = 0.41 


21.08 


18.18 


15.58 


24.22 


20.96 


18.04 


ujb = 0.40 


22.43 


19.33 


16.55 


25.78 


22.32 


19.20 


CUB = 0.38 


25.56 


21.97 


18.77 


29.50 


25.45 


21.84 



2X'^Ri,rsinSsin(j)3 



(l + {X'^Rbrf +2X^RbrcosScos(l)3) 



(12) 



by using Eq. (|TT|). From Eq. (|T2D, we see that the CPA is associated with both weak CP 
and strong phases. The results as a function of angle 03 are plotted in Figure ^. 

In the following, we present the implications of our results. With the limit of SU(3) 
symmetry, it is known that T{Bd -^ K+n-) ^ r(fi, -^ K+R-) and T{Bd -^ K^K^) ^ 
\VtdlVts\^T{Bs -^ K+K-) g. With including SU(3) breaking effects, one has § 



Br (b, -^ K^K^ 



K^^ 



Br B 




:^BsK 



(0)' 



^F^d^(o); 

Vts F^^^(O) 



BriBd^ K^Ti^ 



VtdF^^^iO) 



Br B. 



Or>0 



K^'K 



(13) 

(14) 
and 



We note that Br {Bd — * K^n^) has been observed in the present B factories 0, |^, 
although there is no limit on Br (Bd -^ K^K^j, its estimation has been done by the PQCD 

in Ref. ^. By taking Br {Bd -^ K^n^) ~ 18.5 x 10"^ and Br (Bd -^ K^K^) ~ 1.4 x 10-^ 

we get Br {B, -^ K^K^) ~ 22.68 x 10"^ and Br (b^ -^ K^K^) ~ 26.05 x 10"^. Comparing 
with the values in Table ^, obviously the results are consistent with the those from the 
MPQCD for ub = 0.4, m?^ = 1.7 GeV and 03 ~ 72°. If so, the MPQCD prefers 03 to be 
less than 90°. On the other hand, if neglecting the small difference from EW effects (~ 7% 
difference in our analysis), from Eq. ([TT|), we find that |^ 



Br (Bs -^ K^K^) ^ Br (Bs -^ K'^K^) (l + 2 X"^ Rbr cos 6 cos i 



(15) 



In Eq. ([T5|) , if the interference term associated with cos03cos(5 is negative, it gives Br{Bs -^ 
K^K^)< Br{Bs -^ K^K^). In the MPQCD, it also implies that 03 < 90°. From Eq. (|T3|), 
we also expect that Acp [Bg -^ K'^K^) ^ Acp {Bd — > K'^ii^). Hence, by the measurements 
of Bg — > KK decays, the sign of cos03cos5 can be determined and the predictions of the 
MPQCD can also be verified. 

Finally, we give a brief remark on the relations between Bg — > K^K^ and Bd — > n'^n^ 
decays. It is known that, as shown in Ref. 0, there are close relationships for transition 
form factors between the decays by interchange of d and s quarks, called the U-spin trans- 
formation. Under the U-spin limit, both decays approximately have the same r and 6 as 
defined in Eq. (|TT[) . Moreover, we find that the term associated with cos 5 cos 03 in the 
Bd — ^ TT^TT^ decay has an opposite sign to that in Bg -^ K^K^. This implies that while one 



6 



of both distributions increases with respect to 03, the other one will decrease [^. Hence, by 
precise measurements on the BRs of B^ -^ tt'^tt^ and Bg — ^ K^K^, we also can obtain the 
information of cos 6 cos 03. 

In summary, we have studied the Bg -^ KK decays with the MPQCD approach. We have 
verified that the typical scale in the MPQCD is around 1.7 GeV and the large absorptive 
parts make sin5 ~ —0.45 so that the CP asymmetry could be as large as 15%. Since the 
Tevatron Run II has started a new physics run at \fs = 2 TeV and will collect a data sample 
of 2 fb~^ in the first two years [^. At its initial phase with lOK of Bg -^ KK events, the 
BRs and CP asymmetry predicted by the MPQCD can be tested. 
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Figure Captions 

Figure 1: Form factors of (a) Bd ^ K and (b) Bs ^ K with respect to m^l^. The 
dashed, sohd and dashed-dot hnes correspond to uob =0.38, 0.4 and 0.41 GeV, 
respectively. 

Figure 2: BRs (in units of lO"*^) of (a) B, -^ K^K^ and (b) 5, -^ K+R- with fixing 
the typical scale on WCs. The dashed, dashed-dot, dotted and dense-dot lines 
stand for t=1.5, 1.7, 2.5, 4.8 GeV, with ub = 0.4 GeV and 03 = 72° while the 
solid line is the result of t as a running scale, respectively. 

Figure 3: The CP average BR of Bg — > K^ K~ as a function of 03 with ojb = 0.4 GeV. 
The solid, dashed and dashed-dot lines denote m^^ =1-7, 1.6 and 1.5 GeV, 
respectively. 

Figure 4: The CP asymmetry of Bg — > K^K^ with m% = 1.7 GeV. The dashed, solid 
and dashed-dot lines represent ub =0.41, 0.4 and 0.38 GeV, respectively. 
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Figure 1: Form factors of (a) Bd ^ K and (b) Bg ^ K with respect to m^j^. The dashed, 
sohd and dashed-dot hnes correspond to ujb =0.38, 0.4 and 0.41 GeV, respectively. 
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Figure 2: BRs (in units of 10"^) of (a) B, -^ K'^K^ and (b) B^ -^ K+R- with fixing the 
typical scale on WCs. The dashed, dashed-dot, dotted and dense-dot lines stand for t=1.5, 
1.7, 2.5, 4.8 GeV, with ub = 0.4 GeV and 03 = 72° while the solid line is the result of t as 
a running scale, respectively. 
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Figure 3: The CP average BR of Bg -^ K^K as a function of 03 with ub = 0.4 GeV. The 
sohd, dashed and dashed-dot hues denote rnPj^ =1-7, 1.6 and 1.5 GeV, respectively. 
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Figure 4: The CP asymmetry of Bs — > K^K with m% = 1.7 GeV. The dashed, sohd and 
dashed-dot lines represent ub =0.41, 0.4 and 0.38 GeV, respectively. 
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